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ABSTRACT 

We show that the gaseous halos of collapsed objects introduce a substantial cu- 
mulative opacity to ionizing radiation, even after the smoothly distributed hydrogen 
in the intergalactic medium has been fully reionized. This opacity causes a delay of 
Az > 1 in redshift between the time of the overlap of ionized bubbles in the intergalac- 
tic medium and the lifting of complete Gunn-Peterson Lya absorption. The minihalos 
responsible for this screening effect are not resolved by existing numerical simulations 
of reionization. 
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1. Introduction 

The spectra of the quasars SDSS 1030+0524 at z = 6.28 (Becker et al. 2002) and SDSS 
1044-0125 at z = 5.73 (Djorgovski et al. 2002) feature broad regions of high optical depth, possibly 
implying the existence of a Gunn-Peterson trough of complete resonant Lya absorption by a neutral 
intergalactic medium (Gunn & Peterson 1965). However, the discovery of a Lya emission line from 
a galaxy at z = 6.56 (Hu et al. 2002) implies that the gas in the vicinity of that galaxy is mostly 
ionized and that the neutral fraction is highly non-uniform. The appearance of a Gunn-Peterson 
trough is generally attributed to absorption by a smooth intergalactic medium (hereafter IGM) with 
a significant neutral fraction, and it indicates that observations are approaching the reionization 
era. 

The reionization of hydrogen is expected to involve several distinct stages. The initial, "pre- 
overlap" stage [using the terminology of Gnedin (2000)] consists of individual ionizing sources 
turning on and ionizing their surroundings. The radiation from the first galaxies must make its 
way through the surrounding gas inside the host halo, then through the high-density region which 
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typically surrounds each halo. Once they emerge, the ionization fronts propagate more easily into 
the low-density voids, leaving behind pockets of neutral, high-density gas. During this period the 
IGM is a two-phase medium characterized by highly-ionized regions separated from neutral regions 
by sharp ionization fronts. The main, relatively rapid "overlap" phase of reionization begins when 
neighboring H II regions start to overlap. Whenever two ionized bubbles are joined, each point 
inside their common boundary becomes exposed to ionizing photons from both sources. Therefore, 
the ionizing intensity inside H II regions rises rapidly, allowing those regions to expand into high- 
density gas. This process leads to a state in which the low-density IGM is highly ionized and ionizing 
radiation reaches everywhere except for gas located inside self-shielded, high-density clouds. The 
ionizing intensity continues to grow during the "post-overlap" phase, as an increasing number of 
ionizing sources becomes visible to every point in the IGM. 

The recent observations of strong Lya absorption at z ~ 6 indicate a significant neutral fraction 
in the IGM (e.g., Fan et al. 2002; Lidz et al. 2002). However, they do not necessarily correspond to 
the overlap stage of reionization, since even after the end of overlap, variations in the gas density 
and the ionizing intensity in the IGM may be able to produce such large regions of absorption 
(Barkana 2002). When the voids become transparent to ionizing photons, the naive expectation is 
for the ionizing intensity to increase quickly as multiple ionizing sources begin to reach every point 
in the IGM. However, such a sudden, homogeneous end to the Gunn-Peterson trough may conflict 
with the observations, which imply that along one line of sight, overlap occurs at z > 5.9 (Becker et 
al. 2002), while a relatively large region of strong absorption is found along a second line of sight at 
z = 5.3 (Djorgovski et al. 2002); furthermore, the recent z = 6.56 galaxy may indicate the existence 
of a transparent region of the IGM at this much higher redshift (Hu et al. 2002). Thus, the ionizing 
intensity during post-overlap must have been inhomogeneous and must have increased gradually, 
with both effects possibly caused by the shadows due to the remaining neutral gas. Much of this 
gas may have been located inside the dense halos of collapsed objects. 

In the popular ACDM cosmology, most of the collapsed gaseous halos prior to reionization 
had a virial temperature below 10 4 K, at which cooling due to atomic transitions is heavily sup- 
pressed [Barkana & Loeb (2001), and references therein]. The gaseous low-mass halos (hereafter 
"minihalos" ) thus remain intact until photoionization heating by the cosmic UV background photo- 
evaporates much of this gas back into the IGM during the overlap phase (Barkana & Loeb 1999). 
Haiman, Abel, Sz Madau (2000) noted that these halos may consume a significant number of ioniz- 
ing photons if they are assumed to be completely photo-ionized throughout the time it takes them 
to expand back into the IGM. However, Cen (2001) pointed out that as an ionizing front approaches 
a dense gas halo, it slows down and is eventually halted outside the virial radius, with a shock front 
entering the gas halo instead. Regardless of the precise dynamics, the process of photoevaporation 
takes an amount of time which is set by the characteristic sound speed of ~ 10 km s^ 1 (correspond- 
ing to T ~ 10 4 K). Whether minihalos play a substantial global role during this timespan depends 
on their ability to screen ionizing sources. 

The limited mass resolution of current numerical simulations of reionization (Gnedin 2000; 
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Razoumov et al. 2002) does not allow them to resolve these minihalos and include their cumulative 
screening effect; Haiman, Abel, & Madau (2000) made a crude estimate of the covering factor 
of minihalos and found that it could be significant. In this paper we gauge the importance of 
absorption by minihalos with a detailed semi-analytic calculation. In §2 we describe our model 
for the basic halo properties and the statistics of their screening. In §3 we describe our numerical 
results. Finally, §4 discusses the main conclusions of this work. Throughout the paper we use the 
ACDM cosmological parameters of Q. m = 0.3, VL\ = 0.7, and f2& = 0.05, for the density parameters 
of matter, cosmological constant, and baryons, respectively. We also assume a Hubble constant 
h = 0.7 in units of 100 km s _1 Mpc _1 , and a primordial scale-invariant (n = 1) power spectrum 
with a value a% = 0.8 for the root-mean-square amplitude of mass fluctuations in spheres of radius 
8 h- 1 Mpc. 



2. Theoretical Model 

2.1. Basic Halo Properties 

We begin with a description of our model for the properties of minihalos and of ionizing sources 
in the context of hierarchical galaxy formation in ACDM. The input for our opacity calculation 
involves the distribution of halo masses for halos that contain uncooled gas and for halos which 
host galaxies, as well as knowledge of the density distribution inside each minihalo. 

We assume that the abundance of halos is given by the model of Press & Schechter (1974), with 
the modifications which Sheth & Tormen (1998) and Jenkins et al. (2001) used in order to fit more 
accurately the mass function seen in numerical simulations. Galaxies form in halos in which gas can 
accumulate and cool. At high redshift, gas can cool efficiently in halos down to a virial temperature 
of ~ 10 4 K or a circular velocity of V c ~ 16.5 km s _1 with atomic cooling, which we assume to be 
the dominant cooling mechanism (we consider the effect of H2 cooling in §3). Before reionization the 
IGM is cold and neutral, and these cooling requirements set the minimum mass for halos which can 
host galaxies. During reionization, the ionized IGM is heated to a temperature Tiqm ~ 1-2 x 10 4 
K, and gas infall into halos is suppressed by the increased pressure force. However, this pressure 
suppression is not expected to cause an immediate suppression of the cosmic star formation rate, 
since even after fresh gas infall is halted the gas already in galaxies continues to produce stars, and 
mergers among already-formed gas disks also trigger star formation. We neglect the suppression 
effect in this paper, since we consider only the era of reionization and its immediate aftermath. 

We use the term 'minihalos' to denote halos in which gas accretes but is unable to cool. These 
halos contain virialized dark matter and gas, where the gas does not collapse onto a disk and does 
not participate in star formation. The maximum total mass of a minihalo is 1.6 x 1O 8 M at z = 6, 
corresponding to the atomic cooling threshold (e.g. Barkana & Loeb 2001). This implies that the 
gas in minihalos is almost entirely neutral prior to reionization. The minimum mass is determined 
by the IGM temperature, which prevents gas from collecting inside the smallest halos. The Jeans 
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mass is the minimum mass-scale of a perturbation for which gravity overcomes pressure, in linear 
perturbation theory. The Jeans mass is related only to the evolution of perturbations at a given 
time, while a collapsing halo is the end result of the long evolution history of a perturbation. When 
the Jeans mass varies during this time, the overall suppression of the growth of perturbations 
depends on a time-averaged Jeans mass. Gnedin & Hui (1998) showed that the correct time- 
averaged mass is the filtering mass Mp, which we adopt as the minimum mass of minihalos. At 
z = 6, for example, we find Mp = 4.6 x 1O 4 M in areas that have not yet been reionized and 



If gas collects inside a halo and also cools, it can collapse to high densities and form stars (or 
possibly also a quasar). The ability of stars to form is determined by gas accretion which, in a 
hierarchical model of structure formation, is driven by mergers of dark matter halos. In order to 
determine the lifetime of a typical source, we first determine the time spent by gas (i.e. the gas 
'age') within a given halo using the average rate of mergers which built up the halo. Based on the 
extended Press-Schechter formalism (Lacey & Cole 1993), for a halo of mass M at redshift z, the 
fraction of the halo mass which by some higher redshift Z2 had already accumulated in halos with 
galaxies is 



where D(z) is the linear growth factor at redshift z, S{M) is the variance on mass scale M (defined 
using the linearly-extrapolated power spectrum at z = 0), and M m \ n (z2) is the minimum halo mass 
for hosting a galaxy at Z2 (determined by cooling as discussed above). 

We estimate the total age of gas in the halo as the time since redshift Z2 where Fm(z, z-i) = 0.2, 
so that most (80%) of the gas in the halo has fallen into galaxies since then. Low-mass halos form 
out of gas which has recently cooled for the first time, while high-mass halos form out of gas which 
has already spent previous time inside small galaxies. We emphasize that the age we have defined 
here is not the formation age of the halo itself, but rather it is an estimate for the total period 
during which the gas which is incorporated in the halo has participated in star formation. However, 
the rate of gas infall is not constant, and even within the galaxy itself the gas likely does not form 
stars at a uniform rate. Indeed, galaxies could go through repeated cycles of a star formation 
burst followed by feedback squelching, followed by another cycle of cooling, fragmentation and star 
formation. The details involve complex astrophysics, which are not understood even at low redshift, 
so we account for the general possibility of bursting sources by adding a parameter £ < 1, the duty 
cycle. When £ < 1, compared to £ = 1 there are fewer active sources at any given time (by a factor 
of £) but each has a larger star formation rate (by a factor of 1/C)- I n addition, the star formation 
rate is proportional to the efficiency parameter 77, which is the fraction of gas in galaxies which 
turns into stars. 



heated. 




(1) 



Specifically, the ionizing source in a halo of mass M emits, over its entire lifetime (which at 
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z = 6 is typically ts ~ 1 x 10 8 (£/0.25) yr), the following number of photons per unit solid angle: 

rfo" b u \kpm q ) V 40 A 0.17 J ■ [ > 

In this equation, iV"i on is the overall number of ionizing photons produced per baryon in galactic 
disks; it includes 4000 ionizing photons produced per baryon in stars [for a Scalo (1998) stellar 
IMF and a metallicity equal to 1/20 of the solar value 1 ], times an efficiency r] = 10%, times an 
escape fraction /e S c k = 10% of ionizing photons from the disk itself [see Barkana Sz Loeb (2001) 
and references therein]. Our choice of typical parameters £ = 0.25 and r] = 10% is loosely based 
on numerical simulations and on the low redshift cosmic star formation rate, respectively, but the 
exact parameters valid during reionization can only be determined from future observations. 



2.2. Screening due to Minihalos 

At any given redshift we consider a single halo emitting ionizing photons, surrounded by a 
distribution of mini-halos of various masses at various distances R from the center of the source 
halo. Except where indicated otherwise, we neglect the ionizing radiation from other sources when 
calculating the transmission of the radiation from a source through the minihalos. This assumption 
is justified as long as we consider the early stage of reionization when the ionizing background has 
not yet been established, and as long as we consider a small enough distance that the ionizing 
intensity is dominated by the central ionizing source. 

In order for the radiation to reach a large distance R, it must overcome two obstacles caused 
by the minihalos: (i) it must ionize the neutral gas within the minihalos, and (ii) it must overcome 
the recombination rate of this gas in order to keep it ionized. Our statistical treatment of both 
quantities is similar, so in general we label the quantity we are considering by v. Thus, for example, 
v\ is the sum of the number of hydrogen atoms per unit solid angle encountered by a line of sight 
going out to radius R. This sum is then compared to the total photon emission per solid angle 
dP/dQ of the source, and the photons escape out to radius R along those directions for which 
dP/dQ > v\. Separately we also calculate f2, the sum of the number of recombinations per second 
of hydrogen atoms (assuming the gas starts out ionized) per unit solid angle along a line of sight 
out to radius R. Photons can only escape along those directions for which d 2 P/[dQ dt] > 1/2, where 
the (average) emission rate of photons per solid angle is 

d 2 P _}_dP m 

dndt t s dn ' u 

in terms of the source lifetime ts- 



lr This low metallicity is expected to roughly characterize z ~ 6 galaxies, which form out of an IGM which may 
have been enriched by an even earlier generation of galaxies; see, e.g., Barkana & Loeb (2001) for further discussion. 
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Both v\ and v 2 are determined by the density distribution of gas in the minihalos. Numerical 
simulations of hierarchical halo formation indicate a roughly universal spherically-averaged density 
profile for dark matter in the resulting halos (Navarro, Frenk, & White 1997, hereafter NFW), 
though with considerable scatter among different halos. Such high-resolution simulations have not 
been done with gas. However, in one-dimensional simulations of dark matter and adiabatic gas, 
we find [based on the methods of Thoul & Weinberg (1996) and Haiman, Thoul, & Loeb (1996)] 
similar final density profiles for the gas and dark matter, where both components form a cusp and 
the gas profile is only slightly shallower than the dark matter profile. Note that this result is found 
even though the physics of the two components is rather different, with the dark matter undergoing 
violent relaxation due to shell crossing, while the gas is halted at a virialization shock which moves 
outward. For simplicity, we therefore assume that even in the three-dimensional case the virialized 
gas profile is similar to that of the dark matter, and so we adopt the radial NFW profile for the 
gas density in minihalos, 

where x = r/r v ; r , cn is the concentration parameter, and the virial radius is 

where in this equation and in what follows we assume a sufficiently high redshift for the Einstein-dc 
Sitter limit [VL m (z) k, 1]. The value of the concentration parameter for minihalos at high redshift 
is unknown. We extrapolate from simulations of Bullock et al. (2000), which indicate 

For the population of minihalos we adopt for simplicity a single value of cn, the value expected for 
a minihalo of mass equal to the mass-weighted mean minihalo mass. 

Consider a line of sight which passes through a minihalo of mass M, with an impact parameter 
of a (in units of r v - 1T ) relative to the center of the minihalo. The total hydrogen column density is 
then 

where we include gas out to a maximum radius x m (in units of r v ; r ), and where 



0.7, 

(7) 



5 ( CN ) = log(l + c N )- T ^- • (8) 

1 + CN 

Defining 

CO! = 1 - a 2 c N (9) 

and 

uj 2 = 1 + c N x m , (10) 
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we derive the function 



/i(a,c N ,x m ) = - 
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Similarly, the total recombination rate per unit area, along the same line of sight, is 



r rec = 5.4xl0 4 |V /2(«,c N ,x m ) I 
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Using wi and LO2 from above and defining 
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we find 
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2.3. Statistics of Screening 

The quantities v\ and z/ 2 introduced in the previous subsection are sums over the contributions 
from all minihalos encountered along a line of sight, where a single cross-section of a minihalo at 
proper distance R contributes NmR 2 to v\ and r rec i? 2 to z^ 2 . We wish to calculate the optical depth 
dr for getting a value in the range v\ to v\ + dv\ by intersecting a minihalo along a line of sight. 
This involves an integral along the line of sight of the minihalo abundance times cross-sectional 
area. We first note that the cross-sectional area of a single minihalo, for a given fractional impact 
parameter a, is wa 2 r 2 ir . Thus, the contribution to a given value of v\ is 2ita da r 2 ir ; determining 
the range da that contributes to a given range dv\, we arrive at the total optical depth: 

dr = j R ™ dR f M ™ dn(M) z) 6(fl> ^ z) ^^(^ z) /i(a,c N ,x m ) ^ (15) 

Jo JM min \g^fl{Oi,C N ,X m )\ 

where dn is the halo abundance, b(R, M, z) is the bias of minihalos, and everywhere in the integrand 
a is determined by inverting the equation Nm(M, z, a)R 2 = v\. We get a similar formula for v 2 as 
well. 
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A crucial component of the minihalo abundance is the bias b(R, M, z), which accounts for 
the biased distribution of minihalos near the source halo. We discuss here the calculation of bias 
in some detail. First, we set the bias equal to zero within the virial radius of the source halo, 
since minihalos cannot exist inside the source halo. Outside the virial radius, the bias can be 
divided into two components. One can be called "Lagrangian bias", or "initial condition bias", 
and results from the fact that the initial perturbations in nearby points are correlated. Although 
various analytical approaches to Lagrangian correlations have been suggested, we use the natural 
and general approach of Scannapieco & Barkana (2002), who worked directly within the standard 
theoretical framework of excursion sets (Bond et al. 1991; Lacey &; Cole 1993). Their results yield 
the bias of two halos at a given distance as a function of the two halo masses and the two formation 
redshifts (although in this application we require only the equal redshift case). Most importantly 
for the current application, this approach to bias explicitly accounts for the existence of the source 
halo. For example, there is a strong correlation between nearby rare peaks in the initial density 
field, which tends to give a positive bias among massive halos at high redshift; on the other hand, 
all halos are anti-biased at initial positions that are very close to a given halo A, since much of 
the mass there is likely to belong to the halo A itself when it finally forms. Because of this effect, 
minihalos are anti-biased at initial positions that are close to a source halo, while source halos 
themselves are positively correlated at the same distances (though both are anti-biased at the very 
smallest distances). 

The second part of the bias can be called "Eulerian bias" or "infall bias", and is a result of 
the high-density environment around the massive source halo when it forms, due to infall of mass 
from a large initial region in the direction of the forming halo. If, for example, minihalos contain 
25% of the mass density in an unbiased region, then in the absence of Lagrangian bias, minihalos 
would contain 25% of the mass density also near the source halo; thus, the absolute number density 
of minihalos would increase in proportion to the general mass overdensity near the source halo. In 
general, the Lagrangian bias must be included as an additional factor beyond this. In the spherical 
approximation, the Eulerian bias depends only on the final density profile near the virialized source 
halo; assuming there is no shell crossing (a reasonable approximation outside the virial radius), we 
can map each mass shell in the final density profile to its initial "Lagrangian" position, and we 
then use this initial position in order to calculate the Lagrangian bias factor. There are a number 
of possible approaches to estimate the density profile outside the virial radius of the source halo, 
but approaches based on spherical collapse, for example, do not match the profile from numerical 
simulations near the virial radius. However, all approaches give three-dimensional density profiles 
with power-law slope between R~ 2 and i?~ 3 , and in any case, realistically we expect a large scatter 
in the profiles around halos. For simplicity, here we extrapolate the NFW profile of the source 
halo to larger radii, an approach which also gives a fall-off between R~ 2 and R~ 3 in the relative 
overdensity (i.e., we assume that at large R the density approaches its cosmic mean value). 

In order to calculate the fraction of ionizing photons from the source that escape out to radius 
R, we require the probability distributions of v\ and v-i- However, the above optical depth dr 
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does not directly yield the probability; it only accounts for the possibility of intersecting a single 
minihalo, while there can be multiple intersections along a single line of sight. In addition, the total 
optical depth can be greater than unity. We determine the actual probability distribution with a 
Monte Carlo approach. Specifically, we divide the total optical depth to ~ 100 bins, each of which 
contains a At « 1. In each Monte Carlo trial, we obtain a total v\ (for example) by summing 
over the bins; a given bin makes a contribution with probability Ar, and the contributed value is 
determined by the distribution dr/dvi within the bin. We determine the probability distribution 
using ~ 100,000 trials (for each set of parameters). Note that our calculation assumes that the 
biasing of the minihalo abundance is determined by the correlation with the nearby source halo, 
and we neglect any additional correlations among the minihalos themselves. 

Once we have obtained the probability distributions of v\ and ^2, we can calculate the proba- 
bility distribution of the escape fraction of photons. We often find it convenient to summarize this 
distribution with a single quantity, the mean escape fraction of photons, averaged over the proba- 
bility distribution (which corresponds to averaging over solid angle). Including only the constraint 
of initial ionization (i.e., the escape fraction is 

(2) 

where dp/dv\ denotes the probability distribution of v\. We similarly get / CS c , with v\ and dP/dQ 
replaced by v-i and d 2 P/[dfl dt], respectively. Combining the two constraints exactly would be 
difficult given the fact that they are correlated, since a given minihalo contributes to both v\ and 

V2\ in practice, we estimate the final escape fraction by taking the lower of the two, which is usually 

f (2) 

J esc • 



2.4. Characteristic Distances 

When we consider the escape of photons out to various distances from a source, some values of 
distance are particularly important. We first want to determine a fiducial distance R that represents 
the beginning of the overlap stage of reionization. The overlap stage starts once neighboring H II 
regions begin to merge, thus increasing the ionizing intensity within these regions because of the 
access to multiple ionizing sources. The rising intensity allows the H II regions to expand faster 
and into denser regions, leading to additional mergers, with the process accelerating until the entire 
IGM becomes highly ionized except for gas located inside self-shielded, high-density clouds. 

One guess for the distance R which characterizes the beginning of overlap is simply the typical 
distance to the nearest other ionizing source. However, this value of R does not properly represent 
the beginning of overlap; first, it may be too large since the photons need not reach the nearby 
source itself but only the H II region produced by that source; and second, this R may be too small, 
since if the nearby source has a mass much smaller than that of the first source then it will make 
a negligible addition to the ionizing intensity and this will not advance the process of overlap. We 
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thus choose a distance R which more accurately reflects these considerations. First we find that 
distance where we would typically expect the total mass in additional ionizing sources to be equal 
to the mass M of the first source 2 . Note that we only count active sources (when the duty cycle 
C < 1). Now, since the other sources also send photons toward the first source, we take only a 
fraction of the distance (where just the portion beyond the virial radius is reduced). The precise 
fraction depends on the actual masses of the additional sources; we adopt a fraction of 1/2, which 
is exact for the symmetric case of having a second source halo of mass equal to M. This final choice 
of proper distance R, which represents the beginning of overlap, is referred to hereafter as half the 
equal-mass distance. 

A second important distance is related to the Gunn-Peterson trough. The essential feature 
that defines Gunn-Peterson absorption is the presence of long stretches of continuous, strong 
Lya absorption. Such absorption has been revealed in recent observations at z ~ 6 (Becker et 
al. 2002; Djorgovski et al. 2002). Neutral hydrogen absorbs photons near the Lya resonance with a 
very high opacity as given by Gunn <fc Peterson (1965). As a result, the IGM itself produces strong 
absorption unless it is highly ionized. The IGM is highly ionized only if the ionizing intensity 
is sufficiently high. For a given distribution of ionizing sources, each point in the IGM must see 
sources out to some distance R in order to reach the necessary total ionizing intensity. This proper 
distance R is given by [see equation (18) of Barkana (2002)] 

Rgp = 3A {^) {—) Us J U^J UrwW Uo-otJ UzJ Mpc - 

(17) 

We have defined here the end of the Gunn-Peterson trough as the time when the optical depth 
drops below TL yQ - = 2.5. The total gas fraction inside source halos is F co i. For the purposes 
of this estimate we have also used the approximation that all sources have the same lifetime ts, 
where everywhere in this equation we have assumed £ = 1 since the value of £ does not affect 
the total ionizing intensity. The density in voids which typically dominates the Lya transmission 
is A r = 0.4 in units of the cosmic mean gas density. We have neglected source clustering, but 
such clustering only makes voids relatively empty of sources, increasing the distance R required to 
prevent transmission. For additional discussion of equation (17), see Barkana (2002). 

When we investigate the end of the Gunn-Peterson trough, we must also consider the effect 
of the ionizing background in partially ionizing the minihalos. This is required because the end of 
the Gunn-Peterson trough occurs relatively late, after an ionizing background may have already 
been established; furthermore, it is difficult for the central ionizing source to dominate the ionizing 
intensity at the rather large distance -Rgp- We estimate the effect of the ionizing background as 
follows. First, we find the distance i?o.5 at which / e sc = 0.5 for a typical source (i.e., a source 



2 We avoid the case where the average mass M corresponds to a halo of mass M2 S> M , which is present only rarely 
(with probability M/M2); thus, we choose the distance where we expect to either find a single halo with M2 > M, 
or a total mass M in source halos for which M2 < M individually. 



- 11 - 



with mass equal to the mass- weighted mean mass of all sources). We assume that the ionizing 
background at a random point in the IGM is given by the average integrated intensity due to 
all sources out to the distance R0.5, times one half (since / esc = 0.5 at this distance). With this 
background intensity, ionization equilibrium yields the neutral fraction at each radius of a given 
minihalo. For simplicity, we find the neutral fraction at the virial radius of a typical minihalo, and 
we then assume this neutral fraction to apply throughout every minihalo. This overestimates the 
effect of the ionizing background, since the central portions of minihalos are in reality denser and 
more highly neutral; however, as shown in §3, even with this overestimate the overall effect of the 
partial ionization on screening is small. Once partial ionization is included, / esc increases at a given 
R, and thus i?o.5 also increases; we iterate in order to find the self-consistent value of R0.5, and the 
corresponding value of the neutral fraction in minihalos. We then use this neutral fraction when 
we compute screening at R = Rqp. 

2.5. Direct Gunn-Peterson Absorption 

In addition to their screening of ionizing photons, minihalos also absorb Lya photons and may 
therefore contribute to the Gunn-Peterson trough. Direct absorption near the Lya resonance is 
particularly strong (Gunn & Peterson 1965; see also the related discussion in the previous sub- 
section). However, minihalos contain ~ 20% of the gas in the universe at an average overdensity 
of ~ 200; since for any population of absorbers the covered fraction of a line of sight equals the 
volume fraction of the population 3 , minihalos directly cover only ~ 0.1% of a random line of sight 
to a quasar. More interesting, therefore, is the possibility of absorption due to the damping wings 
of the Lya line (Miralda-Escude 1998). 

To estimate the probability that the population of minihalos produces strong absorption, we 
first consider the optical depth due to a single minihalo; we do this as a function of distance R from 
the minihalo, where the distance represents different wavelengths A by the correspondence of A and 
z via a pure Hubble velocity flow [see the similar approach, e.g., in Loeb & Rybicki (1999)]. We first 
note a number of simplifications: the size of each halo (typically a few kpc) is much smaller than 
the region of significant absorption due to the halo's damping wings (typically a fraction of a Mpc); 
the latter region, in turn, is much smaller than the horizon at the reionization redshift. Under 
these conditions, the damping wing produces an optical depth T"L y a which is simply proportional 
to the total column density of H I divided by the distance squared from the minihalo. This results 
from the Lorentzian damping profile of the line, and it corresponds to the well-known result that 
when the damping wings dominate, the equivalent width of an absorption line is proportional to 



3 More precisely, the covering factor of a line of sight equals the volume filling factor. Consider, for simplicity, a 
number density n of halos which all have the same shape and inclination to the line of sight. A single halo has a 
cross-sectional area A = J A dxdy, where we use two-dimensional (x,y) coordinates. Denote the total depth of the 
halo at each point D(x,y). On a random line of sight, the halos cover a fraction — nj A dx dy D(x,y). This integral 
equals the volume of a single halo. This proof is easily generalized to a population of halos with various shapes. 
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the square root of the column density (note that at small distances from the absorber, distance 
is proportional to change of wavelength AA). Specifically, at a distance R from a column density 

Including the damping wings of the entire population of minihalos, we calculate the total optical 
depth (or covering factor) C along the line of sight, for getting a Lya optical depth above some 
value TL ya . The cross-sectional area contributed by the impact parameter range da is 2iradar^ iT ; 
for a given a, if we determine the maximum distance R from the minihalo for which the optical 
depth is greater than TL ya , then the minihalo covers with this optical depth a total line-of-sight 
distance 2R (including the blue and red damping wings). Thus, the covering factor equals 

ax 

C{>TL ya )= da dn(M,z) 4irar^ ir (M, z) R , (19) 

Jo J M min 

where R in the integrand is determined from T^ ya by inverting equation (18) for each M and a, 
using equation (7) for Nm- As noted in §2.2, we also adopt a single concentration parameter cn 
for all minihalos; this gives them similar density profiles, which simplifies the calculation (i.e., the 
double integral in equation (19) can then be factored into two single- variable integrals). Note that 
the covering factor is inversely proportional to ^Jt Ly a . 

As in §2.3, here too the optical depth includes only the contribution from one minihalo at 
a time, and we determine the full probability distribution of values of TL ya with Monte Carlo 
trials. Note that in this calculation, we do not include clustering among the minihalos themselves. 
Clustering will change the probability distribution somewhat, increasing the optical depth in regions 
of clustered halos while decreasing the optical depth in relatively empty regions. However, the 
damping wings effectively smooth the absorption on the scale of a few comoving Mpc, reducing 
the effect of clustering. There is a similar smoothing of the effect of peculiar velocities, which we 
likewise neglect. 

2.6. Photoevaporation of Minihalos 

As discussed above, the minihalos can affect the evolution of the IGM in a number of ways; 
the gas, however, does not remain in minihalos forever. Barkana & Loeb (1999) showed that 
photoionization heating by the cosmic UV background can evaporate much of this gas back into 
the IGM during the overlap era of reionization. This work also showed that this process affects a 
broad range of halo masses, with only a small gas fraction evaporating out of 10 8 Mq halos, but 
with halos below ~ 1O 6 M losing their entire gas content because of their shallow gravitational 
potential wells. Even after the initial loss of unbound gas, the bound gas in the minihalo expands 
and may further evaporate; however, this stage proceeds more slowly. 
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In this paper we consider two main stages in this process. In the first stage, before reionization 
has begun, the minihalos are fully intact. In the second stage, only the bound gas remains, a 
situation which we approximate by keeping only the minihalos in which > 50% of the gas is 
initially bound after UV heating; this gives a minimum circular velocity V c ~ 11 km s _1 which 
is independent of redshift (Barkana & Loeb 1999). In the transition from the first to the second 
stage, the typical mass of the minihalos that photoevaporate (i.e., their mass- weighted mean) is 
7 x 10 6 Mq. Such a halo has a virial radius of 0.9 kpc, and once it is heated it should evaporate 
with a 10 km s _1 outflow in ~ 90 million years, which is equivalent to the time between z = 6.5 
and z = 6. Additional evolution during the second stage is even slower; the expansion of the 
remaining minihalos (typical M ~ 9 x 10 7 M Q ) takes place on a timescale of ~ 200 million years 
(equivalent to the redshift interval between z = 7.2 and z = 6), and during this time each minihalo 
is expected to accrete additional dark matter or merge with other minihalos. In any case, it is clear 
that if minihalos obstruct reionization then their relatively slow photoevaporation can stretch the 
reionization era over a significant redshift range, as we discuss further below. 

The two-stage process that we have defined is only an approximation and it neglects some 
complications. For example, when an ionizing source begins to ionize the minihalos around it, 
those minihalos can immediately begin to photo-evaporate; however, as just discussed, most of the 
gas in minihalos takes a substantial amount of time to photoevaporate, even in the first stage which 
involves the unbound gas. Our calculations apply during this period. 

3. Results 

Our main goal is to determine quantitatively the degree of effective screening of ionizing sources 
due to the minihalo population. Figure 1 shows the escape fraction / csc of ionizing photons as a 
function of the physical distance R from the source. We pick a representative of the population 
of ionizing sources, namely a halo of mass M = 2.2 x 1O 1O M , which is the mass-weighted mean 
source halo mass (assuming that the minimum halo mass for sources is 1.6 x 10 8 Mq as determined 
by atomic cooling). As shown in the figure, the typical distance from this halo to the nearest other 
ionizing source is R = 24 kpc, giving a high / esc = 92%. However, as discussed in §2.4, the distance 
which best represents the beginning of overlap is half the equal-mass distance. For the case shown 
in the figure, this distance is R = 73 kpc, giving a low / esc = 18%. 

Figure 1 also illustrates the importance of a full calculation of the biased distribution of mini- 
halos and of other source halos around the first source halo. If we did not include bias, we would 
calculate half the equal-mass distance to be R = 0.13 Mpc, and / e sc = 15%. Including only the 
infall bias (see §2.3) increases the abundance of both minihalos and source halos 4 . However, R 



4 When the second halo mass M2 > M, infall bias is computed using infall onto M2. Also, regardless of the choice 
of bias, M is not allowed to fall within the virial radius of M2. 
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Fig. 1. — / esc versus physical distance R from source halo. All curves assume z = 6 and the 
mass-weighted mean source halo mass M = 2.2 x 1O 1O M . Shown is / esc 

assuming no bias (dotted 

curve), infall bias (dashed curve), and full bias (solid curve). On each curve, two distances are 
marked with open circles; these are the distance to the nearest ionizing source (left) and half the 
equal-mass distance (right; see §2.4). 
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is reduced only slightly to 0.12 Mpc, since at this large distance from the virial radius, the total 
enclosed mass is essentially unchanged by infall; however, / csc = 0.03% is still strongly reduced 
at these large distances because of the screening due to the large concentration of minihalos very 
close to the source. To better understand this, consider moving a given minihalo twice closer to 
the source; the resulting halo then contributes 4 times the covering factor, but produces values of 
v\ and i>2 that are reduced by 1/4. We indeed find that infall bias increases the covering factor 
most significantly at the lowest values of u\ and V2 • Finally, when Lagrangian bias is included, R is 
further reduced due to the strong positive correlation among source halos; / esc , however, actually 
increases (at a given R) because of a slight anti-bias of minihalos around the source halo (see §2.3). 

Figure 2 shows the full probability distribution of v>2, from which / csc (shown in Figure 1) is 
derived using equation (16). The figure compares the results of the Monte Carlo calculation, which 
includes the effect of screening due to multiple minihalos, to the naive probability P(> v<i) = 1 — e _T , 
given an optical depth r for having a recombination rate per unit solid angle greater than U2 ■ Clearly, 
the effect of multiple screening changes completely the probability distribution at high values of 
U2] however, whenever the probability of screening the source is significant, it is dominated by 
the contribution of relatively low values of z^, and the mean / csc is only moderately affected by 
multiple scattering. For example, for the three cases shown in the figure, the naive probability 
gives / esc = 5.8%, 21%, and 93% (in order of increasing R), while the exact probability distribution 
yields f esc = 4.2%, 18%, and 92%, respectively. Note also that our inclusion of partial blocking as 
in equation (16) is crucial for getting the correct value of / esc . For example, using the total covering 
factor of all values of U2, the naive probability of having no blocking at all is 2.0%, 7.7%, and 64%, 
respectively; this probability is not an accurate estimate for / esc . 

Although Figure 1 illustrates the results for a single representative source halo mass, a wide 
range of halo masses may in fact contribute to reionization. Figure 3 illustrates how the effective- 
ness of screening varies with halo mass. The variation with halo mass is complicated by several 
competing effects. On the one hand, larger halos produce a higher flux of ionizing photons. On 
the other hand, larger halos must in general send their photons farther before they can effectively 
participate in overlap, i.e., there is an increase with halo mass in the equal-mass distance; however, 
the rate of increase declines at the high-mass end because of the powerful clustering of extremely 
massive halos. The end result is that for the case of z = 6, minihalos screen all source halos up to 
2 x 1O 12 M0 roughly equally, with 13% < / csc < 27%; although halos with even higher masses are 
not strongly screened, most ionizing photons are produced by halos for which / csc is low. 

Figure 3 also shows the separate escape fractions which result from the obstacles (see §2.2) 
caused by initial ionizations {u\) and by recombinations (1^2) ■ In each case, recombinations give the 
stronger constraint, hence we have generally used recombinations to determine escape fractions. 
The figure shows as well the dependence of / csc on redshift, where we make the individual sources 
more efficient at high redshift, so that the total production rate of ionizing photons remains fixed, 
and thus reionization can occur at the higher redshift. At high redshift, the total mass fraction 
in minihalos is smaller, but the biasing around the source is stronger both for minihalos and for 
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Fig. 2. — Cumulative probability P(> v<i) (lower panel) and corresponding P(< v<i) = [1— P(> U2)] 
(upper panel), versus total recombination rate along the line of sight (^2). Shown are the naive 
probability derived from the covering factor (dashed curves), and the actual probability derived 
with a Monte Carlo approach (solid curves). Also indicated is the rate of production of ionizing 
photons by the source (vertical dotted line). The parameters correspond to the solid curve in 
Figure 1, at three different radii. In increasing order, the first two radii are the left and right points 
indicated on the curve in Figure 1, and the third radius is R = 0.1 Mpc. (Increasing R corresponds 
to increasing P(> 1^2) and decreasing P(< ^2)-) 
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Fig. 3. — / csc versus source halo mass M. Shown is / esc at half the equal-mass distance (see §2.4), 
for the populations of sources and minihalos at z = 6 (solid curves), z = 10 (dashed curves), 
and z = 15 (dotted curves). In each pair, the upper curve corresponds to v\ (constraint of initial 
ionization), and the lower curve corresponds to V2 (constraint of recombinations). An open circle 
on each curve marks the mass- weighted mean source halo mass, where the minimum source halo 
mass is assumed to be determined by atomic cooling. All curves assume our case of full bias. 
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other sources. The final result is that for a typical source (whose mass declines with redshift), / esc 
increases with redshift, from 18% at z = 6 to 28% at z = 10 and 44% at z = 15. 

We have thus far adopted the minimum source mass of 1.6 x 1O 8 M given by atomic cooling, 
and have assumed that all halos above this mass produce ionizing photons with the same overall 
efficiency. However, internal stellar feedback effects may be especially disruptive in source halos with 
relatively low masses and low binding energies. Therefore, we consider the possibility of requiring 
a higher minimum mass for halos that host efficient ionizing sources. At the same time we make 
the individual sources more efficient, so that the reionization redshift remains approximately fixed. 
Figure 4 shows that the resulting / esc , for the characteristic source halo mass, increases at a given 
R as the minimum source halo mass is increased. The ionizing photons, however, must reach a 
more distant R before overlap begins, since the abundance of ionizing sources is reduced. The final 
result is that the screening depends only mildly on the minimum source halo mass. For the four 
cases between 1.6 x 1O 8 M and 1.6 x 10 n M o shown in Figure 4, / esc = 18%, 14%, 11%, and 9.2%, 
respectively, i.e., / csc is reduced by only a factor of two when the minimum halo mass is increased 
by three orders of magnitude. 

Clearly, screening due to fully intact minihalos is a substantial barrier and can delay the 
beginning of overlap. Early on, when the dense gas surrounding each minihalo is still neutral, the 
screening is even more effective. For example, when the minimum source mass is given by atomic 
cooling, / csc = 18% if we include the gas out to the virial radius of each minihalo; if we extrapolate 
the NFW profile out to twice the virial radius, around each minihalo, the result is / esc = 4.1%. Thus, 
the surroundings of the minihalos will be ionized before overlap can make substantial progress, and 
at that point the minihalos themselves will begin to photoevaporate. Figure 4 also shows the results 
for the later stage discussed in §2.6, when only the bound gas remains in the minihalos. Once the 
unbound gas has photoevaporated, the effect of screening no longer substantially impedes the start 
of overlap. In particular, the four cases shown in Figure 4 give / csc = 84%, 82%, 78%, and 76%, 
respectively, in order of increasing minimum halo mass. 

The results also depend on additional parameters. For example, our standard case of £ = 0.25 
and rj = 10% gives / e sc = 18% for intact minihalos and 84% for bound minihalos only. Increasing 
C to 0.5 gives 38% and 90% for these two cases, respectively. Note that increasing £ at a fixed 
efficiency makes each source dimmer but increases the number density of sources, thus reducing 
the distance that photons are required to reach in order for the H II bubbles of nearby sources 
to overlap. Reducing ( to 0.25 gives 7.1% and 77%. The results are less sensitive to changing rj; 
setting 7) = 20% gives 24% and 86%, while setting rj = 5% gives 13% and 83%, respectively. Thus, 
the detailed numbers depend on the uncertain source parameters, but the qualitative conclusion 
remains the same: intact minihalos substantially impede the overlap stage, while bound minihalos 
do not. 

Although overlap can proceed once only bound minihalos remain, observationally the universe 
remains opaque to Lya radiation. In other words, the ionizing intensity in the IGM is initially 
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Fig. 4. — / esc versus physical distance R from source halo. Shown are / esc for all mini-halos (solid 
curves), and for mini-halos in which most of the gas remains bound (dashed curves; see text); 
also indicated is half of the equal- mass distance (dotted lines). The four cases assume a minimum 
source mass of 1.6 x 1O 8 M , 1.6 x 1O 9 M , 1.6 x 1O 1O M , and 1.6 x 1O U M , respectively from left 
to right. Each curve is calculated for a halo of mass equal to the mass- weighted mean source halo 
mass (respective values 2.2 x 1O 1O M , 4.1 x 1O 1O M , 1.1 x lO n M , and 4.1 x lO n M ). All curves 
assume z = 6 and our case of full bias. 
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low, and the small remaining H I fraction continues to produce a Gunn-Peterson trough. The 
trough disappears only when sources can send their photons out to the rather large distance Rqp 
(see §2.4). Even the bound minihalos on their own produce (for our standard parameters, and a 
typical source) a covering factor of 71 for total blocking, which implies absolutely total screening 
for R = Rqp- Note that we include in this calculation the effect of the ionizing background in 
partially ionizing the minihalos (see §2.4). The effect of this on the overall screening is small since 
in the current context, when a line of sight from a source encounters a minihalo the typical result is 
complete blocking even though the neutral fraction at the virial radius of the minihalo is ~ 1%; this 
results from the fact that we are considering only bound minihalos (which are relatively massive, 
and thus have high column densities), and we are considering large distances from the source (so 
that the column density or recombination rate per unit solid angle is high). 

Thus, the screening effect can substantially delay the end of the Gunn-Peterson trough, and 
this remains true even if the gas in some minihalos cools due to H2 molecules and only a small 
fraction of the minihalos are left intact. The history of molecular cooling in minihalos prior to 
reionization involves complicated effects of radiative feedback, effects which are now beginning to be 
explored with a new generation of numerical simulation codes. Although only redshifts > 10 can be 
computed and the results must be considered tentative, Ricotti, Gnedin, & Shull (2002) found that 
molecular cooling does occur in some minihalos, but the total ionizing intensity is limited by a strong 
negative feedback of direct ionization of H2 molecules; in particular, they argued that minihalos 
cannot contribute substantially to overlap compared to more massive ionizing sources. Future 
simulations may determine the precise fraction of minihalos in which the gas remains uncooled. We 
find that even if only 5% of the minihalos are intact, the bound ones produce essentially complete 
blocking (i.e., / esc = 2.5%) at R = Rqp] with 2% of the minihalos the escape fraction is still small 
(23%), while 1% of the bound minihalos block around half the photons (i.e., / esc = 48%). 

Finally, we consider the absorption of Lya photons by minihalos (§2.5). Since the covering 
factor (or optical depth) C is a power law in TL ya , the conversion via Monte Carlo to a cumulative 
probability P(> TL yQ! ) involves a universal function, where P is a function of C only and does not 
depend separately on TL ya . This universal function is plotted in Figure 5, where this function is 
shown to be significantly different from the naive probability of 1 — e~ c . The probability is 50% 
when C = 0.53, 68% when C = 0.78, 95% when C = 1.5, and 99% when C = 2.0. 

For the covering factor itself we obtain with fully intact minihalos, C(> r^ ya ) = 0.97/^/TL yQ . 
Therefore, P(> 1) = 78% of the line of sight is covered with TL ya > 1 by minihalo damping wings; 
the number is P(> 5) = 42% with T^ ya > 5. In the earliest stage, when we can extrapolate 
the NFW profile out to twice the virial radius around each minihalo, the corresponding results 
are C(> r Lya ) = with P(> 1) = 99.6% and P(> 5) = 81%. On the other hand, 

in the late stage when only the bound gas remains inside minihalos, C(> TL ya ) = 0.083/^/TL yQ! , 
with P(> 1) = 8.3% and P(> 5) = 3.7%. The low optical depth produced by the bound gas 
results from the fact that, for a given total mass, small halos are more effective Lya absorbers 
than large halos. This results from equation (18): a reduction in column density by a factor of 
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Fig. 5. — Cumulative probability P(> 7L ya ) for Lya absorption with optical depth above TL yQ , 
versus covering factor C(> TL ya ). Shown are the naive probability derived from the covering factor 
(dashed curve), and the actual probability derived with a Monte Carlo approach (solid curve). 
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2 reduces the distance covered (by a given TL yQ ,) by only a factor of \[2. We conclude that as 
long as the minihalos are intact, they can produce significant Lya absorption, although they are 
unlikely to themselves cause full absorption (i.e., TL yQ > 5 with a probability near unity). After 
the unbound gas photoevaporates, however, the minihalos can make only a minor contribution to 
Lya absorption. 

4. Conclusions 

In this paper we have shown that the population of gaseous minihalos can produce substantial 
screening of ionizing photons during reionization. After accounting for the biased distribution of 
minihalos around sources, we have found that a fully intact population can delay the start of 
overlap; for a typical source at redshift z = 6, advancing overlap requires getting the ionizing 
photons out to a proper distance R = 73 kpc, for which the minihalos limit the escape fraction to 
a low value of / esc = 18%. The escape fraction out to the appropriate distance is roughly the same 
for source halos over a wide range of mass, with only rare M > 1O 13 M halos having a much larger 
escape fraction. Once the smallest minihalos photo-evaporate, gas remains in the bound minihalos 
for a substantially longer period. The bound minihalos are unable to delay the start of overlap, 
producing a value of / csc = 84% (for our standard parameters at z = 6). 

To obtain the escape fraction, we have included the effect of screening due to multiple minihalos 
(using a Monte Carlo approach), and also the effect of partial blocking [see equation (16)]. We 
typically adopted the minimum source mass of 1.6 x 10 8 Mq given by atomic cooling; however, if 
stellar feedback favors higher-mass halos, we have shown that / esc is insensitive to the minimum 
source halo mass if the reionization redshift is held fixed; and although the detailed numbers depend 
on uncertain source parameters (such as efficiency and duty cycle) , we have shown that the overall 
conclusion remains the same: the full population of minihalos can substantially impede the overlap 
stage, while the bound minihalos themselves cannot. 

Even after overlap, screening due to minihalos remains highly effective until the minihalos 
photo-evaporate; this screening keeps the ionizing intensity in the IGM low, and the Gunn-Peterson 
trough remains until even the bound minihalos are disrupted. This can produce a lag of ~ 200 
Myr (equivalent to redshift interval z = 7.2-6) between overlap and the end of Gunn-Peterson 
absorption. This screening is so strong that it is highly effective even if molecular cooling occurs 
in most minihalos and only a few percent of them retain uncooled gas. 

Intact minihalos may also directly absorb Lya photons. For the full population at z = 6, 42% 
of a random line of sight is covered with r^ ya > 5. However, as noted above, a full population of 
minihalos screens ionizing sources and makes the IGM absorb Lya photons much stronger than do 
the minihalos themselves. When only the bound minihalos remain, the ionizing intensity is still 
low in the IGM, but at this stage the direct absorption due to minihalos is already insignificant 
(e.g., only 3.7% of the line of sight is covered at TL ya > 5). 
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Current numerical simulations cannot resolve the population of minihalos during the reioniza- 
tion era. However, we have demonstrated in this paper that any conclusions regarding the timing 
of the overlap era of reionization and of the lifting of Gunn- Peterson absorption depend crucially 
on the behavior of the minihalos. 

This work was supported in part by NASA grants NAG 5-7039, 5-7768, and NSF grants AST- 
9900877, AST-0071019 for AL. RB acknowledges the support of an Alon Fellowship at Tel Aviv 
University. 
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